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ABSTRACT 

In the applications of solar magneto-seismology (SMS), employing the ratio of the period of the fundamental mode 
to twice the one of its first overtone, P1/2P2, plays an important role. We examine how field-aligned flows affect 
the dispersion properties, and hence the period ratios, of standing modes supported by magnetic slabs in the solar 
atmosphere. We numerically solve the dispersion relations and devise a graphic means to construct standing modes. 
For coronal slabs, we find that the flow effects are significant, for the fast kink and sausage modes alike. For the kink 
ones, they may reduce P1/2P2 by up to 23% compared with the static case, and the minimum allowed P1/2P2 can fall 
below the lower limit analytically derived for static slabs. For the sausage modes, while introducing the flow reduces 
P1/2P2 by typically < 5% relative to the static case, it significantly increases the threshold aspect ratio only above 
which standing sausage modes can be supported, meaning that their detectability is restricted to even wider slabs. In 
the case of photospheric slabs, the flow effect is not as strong. However, standing modes are distinct from the coronal 
case in that standing kink modes show a P1/2P2 that deviates from unity even for a zero- width slab, while standing 
sausage modes no longer suffer from a threshold aspect ratio. We conclude that transverse structuring in plasma 
density and flow speed should be considered in seismological applications of multiple periodicities to solar atmospheric 
structures. 

Subject headings: magnetohydrodynamics (MHD) - Sun: corona - Sun: oscillations 
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1. INTRODUCTION 



A 



rich variety of magnetohydrodynamic (MHD) 
waves and oscillations have been detected in the 
structured solar atmosphere (e.g., the 

, r-: 5 ; ■ ' ■» tt-^ — — 

by 



reviews 



Aschwandcn 2004'; Nakariakov & Verwichte '2005 
Nakariakov & Erdclyi 2009; Erdclvi & Goosscns 2011 
De Moortel fc Nakariakovl I2012i) . While identifying 



waves/oscillations was possible with radio and op tical 
instruments prior to the 1990s ([Roberts et al.l[T984 and 
references therein), the majority of such observations 
emerged only after the advent of the Solar and He- 
liospheric Observatory (SOHO), the Transition Region 
and Corona Explorer (TRACE), and later Hinode, in 
conjunction with the availability of ground-based in- 
struments with very high temporal and spatial reso- 
lution such as the Dutch Open Telescope (DOT) and 
the Swedish Solar Telescope (SST). For instance, imag- 
ing and spectroscopic instruments on SOHO, along with 
the Hq, measurements with SST and DOT have shown 
that MHD oscillations abound in promi nences - cool 
and d ense clouds suspended in the corona ()Arregui et al.l 
120121 and references therein). Instruments with high spa- 
tial resolution like the Solar Optical Telescope (SOT) on 
board Hinode have detected transverse oscillatory loop- 
like structur es accompanied by downflow s termed "coro- 
nal rains" (jAntolin fc Verwichtd [20TTI ). Hinode/SOT 
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has also detected transverse oscillations in spicules - 
chromospheric materials isolated from a coron al envi- 
ronment [De P onticu ct al. 2007: iHe et"alll2009[ ). Using 
the Ultraviole t Coronagraph Spectrometer (UVCS) data, 
[Morgan et al.l |2004i) detected significant periodic oscil- 
lations in H I Lya intensities in a coronal hole, the quiet 
Sun and a streamer. Coronal loops offer by far the most 
samples of MHD oscillations, among which the major- 
ity seems to be standing transverse oscillations tri ggered 
by fla res (e.g.. .Aschwanden et al..,1999; .Nakariakov et al.l 
IT99I . 

Identifying the observed oscillations with a particu- 
lar MHD mode collectively supported by a magnetized 
structure is physically crucial but proves difficult. This 
is best illustrated by the substantial controversy asso- 
ciated with interpreting the observed transverse oscil- 
lations in t erms of Alfven waves (for the most recent 
review, see IMathioudakis et al.l 1201^. T his happens 
with both coronal (Tomczvk et al. 1 20071 ) and spicule 
measurements (De Ponticu ct al. 200'^), the interpreta- 
tion of which as Alfven waves was suggeste d to b e in- 
consiste nt with theory by Erdelyi fc FedunI (|2007D and 
later by I Van Doorsselaere et al.l ()2008al ). The first un- 
ambiguous detect ion of Alfv e n way es in the solar atmo- 
sphere is due to iJess et all ()2009D where Ha observa- 
tions with SST of a group of magnetic bright points are 
analyzed. Less controversial is the identification of fast 
and slow waves. Consider first the coronal case. Flare- 
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triggered transverse oscillation s supported by loops are 
usually fast kink mode s Ce.g.. lAschwanden et alj 11999 : 
Nakariakov et al.' '19991; iVan Doorsselaere et all 12008b : 
Verwichtc et al. 2009). On the other hand, fast 
sausage mode s were also inferred alb eit much less abun- 
dant (Nakaria kov et al.li2003: Melnikov et al.l 1200 5). As 
for the slow waves, they are frequ ently observ ed and ap- 
pear both in the form of standing ()Wans!jl20Tll and refer- 
ences therein) and as propagating ones (see iDe Moortell 
[200l iBanerjee et aP 120071 : iDe Moortei 120091 for com- 
prehensive reviews). When it comes to atmospheric 
parts other than the corona, s ausage waves a.re dir ectly 
observed in th e photosphere (iMorton et al.l 120111) and 
chromosphere ()Morton et al.l I2012D . Note that in the 
fine structures examined in the latter study, these fast 
sausage waves, manifested as intensity oscillations, ap- 
pear simultaneously with fast kink waves, manifested in 
the form of transverse displacements. That periodic lon- 
gitudinal and transverse oscillatory motions may accom- 
pany each other was first shown by lErdelvi fc TarovanI 
()2008[) for Hinode loops. At this point, it should be noted 
that adopting a helioseismic approach by working with 
the power spectra of, say, the Doppler shift time series, 
can greatly improve t he capability fo r wave identification, 
as was advocated bv lTarovan et al.l ()2007[ ). 

Once identified, the measured properties of waves and 
oscillations can then be combined with MHD theory to 
yield parameters of the solar atmosphere that are oth- 
erwise difficult to measure. The idea behind this was 
put forward as early as i n the 1970 's following the pio- 
neerin g work by lUchidal (|1970[ ) and iZaitsev fc Stepanovl 
(I1975D. Along with othe r s, Ro b erts and co-workers (e.g., 
lEdwin fc RobCTtg [1981 [l98l [Roberts et all [1981 ex- 
pressed the idea in the mathematical form now in rou- 
tine use. This powerful diagnostic capability, ini- 
tially termed "coronal seismology" , may be properly 
called "solar magneto-seismology" (SMS), given that 
both observations and seismological applications ex- 
ten d w ell beyond the solar corona (for recent reviews, 
see [Rob erts 2000, 2008; Nakariakov & Verwichtc 2001; 
iRuderman fc Erdelvil 120091) . The parameter that is 
perhaps most frequently inferred is the magnetic field 
strength. For instance, measuring the periodic fiare- 
generated transverse displacements in several TRACE 
loops leads to concr ete estimates of the corona l mag- 
netic field strength ([Nakariakov k OfmanI |2001[ ). This 
practice was also applied to loops o bserved by Hin- 
ode/ S OT (e.g., Erdclvi fc Tarovan 2(]0l lOfman fc Wang! 
[20081: iVan Doorsselaere et al., ,2008b ). and recently by 
the Atmospheric Imaging Assembly (AIA) on board the 
Solar Dynamic Observatory (SDO) ([White fc Verwichtd 
[20T1 . Certainly SMS can offer more. As a matter of 
fact, slow waves running in Active Region loops mea- 
sured by the Solar Terrestrial R elations Observatory 
(STEREO) enabled [Marsh et all ([2009) to infer loop 
temperatures via the measured propagation speeds. Slow 
waves measured by Hinode also helped derive the e ffec- 
tive adiabatic index (|Van Doorsselaere et al.l[2011a[ ). 

Multiple periodicities detected in some oscillating 
structures are playing an increasi ngly important role 
in solar magneto-seismolo gy (see lAndries et al.l 120091 : 
IRuderman fc Erdelvil [20091 for recent reviews) . They 
may have disparate timescales, differing by nearly a fac- 
tor of 10 as shown in the light curve data obtained by 



the LYRA irradiance experiment , indicative of the coex- 
istence of slow and fast modes ([Van Doorsselaere et al.l 
I2011bf) . This coexistence can help derive the plasma 
beta and density contrast of the flare loops where 
the oscillations are found. They may have compa- 
rable timescales, favoring an interpretation of the co- 
existence of a fundamental mode with its overtones. 
Two (e.g.. lVerwichte et al]|2CI04l: [Van Doorsselaere et al.l 
20071) and three periodicities (IDe Moortel fc BradYll2007 



av\ 

van Doorsselaere et al.l [20091 : llnglis fc: Nakariakovl 12009 ) 
have been found. One thing peculiar with the multiple 
periodicities is that the ratio between the period of the 
fundamental and twice the period of its flrst overtone, 
P1/2P2, deviates in general from unity. For instance, 
the loops labeled C and D in the TRACE 17lA images 
on 15 Apri l 2001 exhibit a Pi /2fi of 0.91 and 0.82, re- 
spectively ([Verwichtc et al.l 120041 ) . A further study us- 
ing similar T RACE 17lA loops on 23 Nov 1998 yields a 
P1/2P2 of 0.9 ([Van Doorsselaere et al.ll2007[ ). This is not 
expected for kink modes supported by a longitudinally 
uniform loop with tiny aspect ratios. Hence the deviation 
of P1/2P2 fro m 1 may signify some longitudinal structur- 
ing. Actually lAndries et al.l ([2005D come up with a one- 
to-one look-up diagram where the density scale height 
along the loop in units of the loop length can be read- 
ily deduced once 1 — P1/2P2 is known. Fundamental or 
global sausage modes with their first overtones were also 
found. For example, fiare-associated quasi-periodic pul- 
sations measured with the Nobeyama Radioheliograph 
(NoRH) on 12 Jan 2000 contain a global sausage mode 
with period Pi ~ 14 — 17 s, and a highe r harmonic with 
P2 ^ 8 - 11 s ([Nakariakov et all I2003I : iMelnikov efall 
120051 ). Moreover, high spatial resolution H^ images of 
oscillating cool post-flare loops measured with the So- 
lar Tower Telescope at Aryabhatta Research Institute of 
Observational Sciences (ARIES) yie ld a Pi w 587 s and 
a P2 « 349 s ([Srivastava et al.ll20'08D . 

Flows are ubiqu itous in the solar atmosphere (e.g., 
lAschwandenl 120041) and have been found in oscillating: 
struc tures (e.g.. lOfman fc Wangl 120081 : iSrivastava et al.l 
[2008( 1. Take the corona for instance. The frequently 
found flow speeds, usually < 100 km and hence well 
below the Alfven speed, have only weak effects on the 
periods of the fundamental ki nk mode and its overtone 
supported by thin loops (e.g.. lRuderman|[201Cl[) . There- 
fore it is suggested that the effect of flow may be ne- 
glected in the seismological applications of 1 — P1/2P2 to 
the corona. However, the flow speeds may not be always 
small. Although less frequent, flow speeds reaching the 
Alfvenic range (^ 10'^ km s~^) have been seen associated 
with explosive events ('e.g.. llnnes et al.|[2003l : iHarra et al.l 
l2005f ). Taking into account a siphon flow of the order 
of the Alfven speed may lead to significant revisions to 
the loop parameters deduced from seismology. For the 
standing kink mo de measured with TR ACE and SoHO 
on 15 Sept 2001 ([Verwichte et al.ll2010f ). iTerradas et all 
([201 It ) show that neglecting the flows leads to an under- 
estimation of the loop magnetic field by up to a factor of 
three! Such a revision is important in its own right, and 
also important in the context of coronal heating because 
the energy flux density carried by the transverse waves 
is proportional to the magnetic field strength. 

Given the importance of multiple periodicities in the 
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context of solar magneto-seismology, and the extent to 
which a flow may affect the seismological applications, 
it seems natural to examine further how the flows af- 
fect multiple periodicities from the theoretical perspec- 
tive. In the c ase of cyli ndrical geometry, this has been 
done by RudermanI (|2010 ) in the thin-tube approxima- 
tion. Here we consider a slab geometry first, and leave 
the study of coronal cylinders with finite aspect ratio to 
a future publication. By doing so, we will complement 
and expand the only study concerning the period ratios 
of a coronal slab (Macnamara & Roberts 2011), which 
shows that the transverse structuring may contribute to 
the deviation of the period ratio from 1. This is true 
even for relatively thin slabs, and is more pronounced for 
hig h density contrasts. Ex t endin g the static slab study 
by iMacnamara fc RobertsI ()2011[ ) to incorporate flows, 
we will examine how the flows affect the dispersion prop- 
erties and hence the period ratios of both standing kink 
and sausage modes. 

This paper is organized as follows. In section [2l we 
briefly describe the slab dispersion relation. Section [3] 
deals with coronal slabs by first giving an overview of 
their dispersion diagrams, then describing a graphical 
means to derive the period ratios, and also examining 
in detail how the fiow affects the period ratios for the 
standing kink and sausage modes. The same approach 
is also applied to isolated photospheric slabs, the results 
thus derived are given in section |4l Finally, section [5] 
summarizes the results, ending with some concluding re- 
marks. 

2. SLAB DISPERSION RELATION 

Consider a slab of half-width d with time-independent 
siphon flows. As illustrated in FigdJ the slab is infinite 
in the y— and z— directions, but is bordered by two 
interfaces x = ±d. Let the subscripts and e denote the 
otherwise uniform parameters inside and external to the 
slab, respectively. The background magnetic fields (Bq 
and Be), together with the flow velocities (Uq and Ue), 
are in the z— direction. Let p and p denote the mass 
density and thermal pressure. From the force balance 
condition across the interfaces it follows that 
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Po 



2cl 



2cl + 7^;2 ^ 



(1) 



where 7 = 5/3 is the adiabatic index, c 
adiabatic sound and va 



y/jp/p is the 



^/B^JAttp the Alfven speed. 
To proceed, it proves convenient to introduce the tube 
speed 



^2 



(2) 



where i — 0,e. 

The dispersion relation (DR) for the waves sup- 
ported by a slab with flow has been examined by a 
number of au thors ( e .g., Nakariakoy fc RobcrtSi 119951 : 
Epardcr & Nar avananI 120001: iMiteva et al.l l2003t ) (also 
see Terra-Homem et al. I I2OO3I for the waves supported by 
cylinders with flows). Restricting oneself to propagation 
in the x — z plane, and to linear perturbations, one may 
adopt the ansatz that any perturbation df{x,z;t) to the 



equilibrium f{x) is in the form 

6f{x, z; t) = Re |/(a;) exp [— i [ujt 



kz 



(3) 



where Re(- • • ) means taking the real part of the function. 
The phase speed Wph is defined as Wph = u /k. Of critical 
importance is the variable 



cf - (wph - Uif 










- {vph - Uif 





, (4) 



where i = 0, e. To ensure the waves external to the slab 
are evanescent, one has to require a positive mg, which 
translates into 

d,e < (Vph - Uef < 4,e Or (Wph - Uef < (5) 

where Cni,e — nrin(ce, f^e) and cu.e — niax(ce, WAe)- On 
the other hand, the signs of mg determine the spatial 
profiles of the perturbations in the transverse (i.e., x—) 
direction. For itiq < (mg > 0), the modes are body 
(surface) waves, corresponding to an oscillatory (a spa- 
tially decaying) x— dependence inside the slab. Evi- 
dently, body waves correspond to 

{vph - Uof > c^^g or 4g < (t;ph - Uof < c^^^q, (6) 

whereas surface waves take up the rest of the possible 

Vph- 

Plugging Eq.Q into the linearized, ideal MHD equa- 
tions, one may derive the DR by requiring the following 
boundary conditions be satisfied: First, the transverse 
(i.e., X—) component of the Lagrangian displacement 
is continuous at x = id. Second, so is the total pressure 
Spt- The DR reads | Nakariakov & Roberts 1995) (see 

tiiir 



also lEdwin fc Robert ^ [19821 for its static counterpart) 
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for surface waves, and 
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for body waves (n^ = ^TOq > 0). Furthermore, the upper 
(lower) case is for kink (sausage) modes. Note that the 
absolute value of rrie appears to ensure the waves are 
trapped, i.e., evanescent outside the slab. 

It proves necessary to examine also the relative impor- 
tance of transverse displacement versus the density fiuc- 

tuation. This is evaluated hy X = {p / po) / {^x / d) , 

x—d 

which reads 



X 



d?Tig(aj — kUof 



PT 



[k^cl-{uj~kUoY]dpT/dx 



(9) 



x — d 



For body waves, {dpT/dx)/pT inside the slab equals 
nocot(noa;) for a kink wave, and — notan(noa;) for a 
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sausage one, leading to 



(nod) tan{n()d) 
kink, 



sausage. 

For surface waves, {dpT/dx)/pT inside the slab is 
mo coth(rno2;) for a kink wave, and toq tanh(TOo2;) for a 
sausage one, leading to 

{(mod) tanh(mo(i) 
kink, . 
(mod) coth(mod) ^ > 
sausage. 



X = 



Before presenting the solutions to the dispersion rela- 
tions ^ and ([U, let us demonstrate the following sym- 
metry properties. 

1. If [uph, fc; C/o, t/e] represents a solution to the DR, 
then so does [wph, —k; Uo, Ue]. This is true because 
both uiq and m^ depend on k only via fc^, and 
me appears in the absolute value bars. Therefore 
changing the sign of k changes the sign of mo , and 
consequently simultaneously changes the signs of 
both sides of the DRs. The end result is that, the 
dispersion diagram expressing Vph as a function of 
k is symmetric about the Wph— axis and one needs 
only to consider the 1st and 4th quadrants, as will 
be done in Figs 12] and H] For more details, please 
see appendix. 

2. When the external medium is at rest, if [wph, k; Uq] 
is a solution, then so is [—Vph,k;—Uo] (see Eq.dH) 
with Ue ~ 0). It follows that the dispersion di- 
agram for some positive Uq is a reflection about 
the k— axis of the diagram for —Uq. Actually 
thi s tendency was clea r already in Figs. 4 and 5 
in iTerra-Homem et al.l ()2003D . even though the 
negative flow speed adopted in constructing Fig. 4 
therein is slightly smaller in magnitude than the 
positive one used in generating Fig. 5. More discus- 
sions on this are given in the appendix. 

We stress that if Ue is not zero, then symmetry 2 does 
not hold. The dispersion behavior, or rather the wave 
behavior in general, critically dep ends on the signs of Up 
and Uq (see the recent rev i ew by iTar ovan fc RudermanI 

m(also lAndries et all (|200Cl[)nAndrics fc Goossens 
). 

3. PERIOD RATIOS FOR STANDING MODES SUPPORTED 
BY CORONAL SLABS 

3.1. Overview of Coronal Slab Dispersion Diagrams 

Let us consider first the coronal case, by which we mean 
the ordering VAe > vao > cq > Ce holds. To be specific, 
we choose cq = 1, vao — 4 and Cg = 0.72. Taking the 
temperature pertaining to TRACE 171A(-- 0.96 MK) 
as representative, one has cq « 130 km s^^, and hence 
VAO = 520 km s~^, and Ce ~ 9 3 km s~^, which are 
very c lose to the values derived by iNakariakov fc OfmanI 
(|2001f) . The reference value for the external Alfven speed 
VAe is taken to be 2vao, corresponding to a density con- 
trast po/pe = 3.76 (see Eq.([l]) which shows that when 



v\, po / Pe ~ {vAe/vAo)'^)- Wc will also examine the 
cases where VAe/vAo is 3 and 4, which correspond to a 
Po/pe of 8.41 and 14.9, respectively. 

Figure [2] presents the phase speeds Vph as a function 
of longitudinal wavenumber k for a series oi Uq. From 
top to bottom, FigsllK to [2li correspond to an Mq of 
0,0.8,1.2 and 3.2, respectively. Here Mq = Uq/cq ex- 
presses the internal flow speed in units of the internal 
sound speed. The dashed (solid) curves are for kink 
(sausage) modes. As shown in FiglJlD, the modes are fur- 
ther labeled by combinations of letters b/f+F/S+K/S, 
representing backward or forward, Fast or Slow, Kink 
or Sausage. While "Fast" or "Slow" is determined by 
the magnitude of the phase speeds, "Backward" or "For- 
ward" is determined by the sign of the phase speeds when 
the flow is absent. The number appended to the letters 
denotes the order of occurrence. As such, bFKl repre- 
sents the first branch of backward Fast Kink mode, and 
likewise, fFS2 the second branch of forward Fast Sausage 
mode. On the left of each panel, the characteristic speeds 
external to the slab are given. Although nearly indistin- 
guishable, there are two windows where waves are pro- 
hibited if their phase speeds are such that c^^ < v"^^ < Cg . 
On the right of each panel, characteristic speeds interior 
to the slab are given. When Uph — Uq lies in the inter- 
vals {-VAO, -Co), (-CTO, cto), or (co, vao), the waves are 
surfaces waves. Otherwise, they are body ones. As was 
pointed out by Nakariakov & Roberts (1995) (hereafter 
NR95), and also evident from Figl^l no surface waves are 
introduced by the flow in the coronal case. The rest of the 
section complements and expands the study by NR95. 

Figure [2] shows that the dispersion diagrams demon- 
strate a clear dependence on the flow speed. With in- 
creasing Uo, there is increasingly less symmetry between 
the lower and upper half of the diagrams, with a perfect 
symmetry taking place only for a static slab (Figl2^). 
For the ease of description, let us consider first the slow 
modes, namely the waves with Vp^ of the order of the 
sound speeds. The most obvious variation with varying 
Uq is that the slow modes are shifted upwards, encom- 
passing always the windows (— Co -|- Uo, —cto + Uo) and 
(cto + Uo, Co + Uo)- This is understandable given that for 
kd^O, 



v^^^Uo±CTo\l + -j;^k^d^, (12) 



where 



^_f(j-l/2)7r kink 



sausage 



(13) 



with j = 1, 2, • • • . In the opposite extreme fed 3> 1, one 
may find 



VpY.-Uo±co\l-^^^{kd)-^, (14) 

V ^AO ^ ^0 



where 



_ ]Ti kink 
^ I ( j — 1/2) TT sausage 



(15) 



with j = 1,2,---. In Eqs.([T2]) and ([141), the plus and 
minus signs correspond to the upper and lower band, re- 
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spectively. They also help explain the occurrence of the 
backward slow waves, labeled BSK and BSS, in the up- 
per half of FigsHt and Hi: even though BSK and BSS 
have negative phase speeds in a frame comoving with 
the internal flow, in the current frame they have posi- 
tive speeds as long as Uq/cto > 1 or Mq > 0.97. In 
view of Egs. p^ and p^ . there actually exists an in- 
finite number of branches of slow modes. However, in 
Fig [2] it is neither easy to differentiate these branches, 
or to distinguish between the kink and sausage modes, 
since cq differs only slightly from cto- Given that the 
slow modes are nearly nondispersive, and that the devia- 
tion from unity of the period ratio P1/2P2 is due entirely 
to wave dispersion in the current study, we will not ex- 
amine P1/2P2 associated with slow modes. 

The fast modes are more interesting since they show 
significant dispersion. With increasing Uq, while the 
modes in the upper half of the Wph — k plane become less 
dispersive, the opposite is true for the ones in the lower- 
half: the phase speeds of fast modes are bounded in the 
former by Uq + vao and VAe, and bounded by Uq — vao 
and —VAe in the latter. It can be readily shown that this 
behavior results from the fact that for kd^ 1. 



Vph f^Uo± VAo\ 1 + 



r(M)- 



(16) 



On the other hand, except for the first kink branch, both 
kink and sausage modes have cutoffs, which are given by 



(17) 



where 



A± = 



(Co+^'Ao)[("AeTC^o)2-4o] 



AOl 



For the first kink branch, an analytic approximation to 
Uph in the slender slab limit fed <C 1 can be found, which 
reads 



V%^±VAeV^-V^ikd)^ 



where 



[{vAeTUu?-vl,?{vle-cl) 



(18) 



(19) 



with Wpj^ and v^^^ representing the upper and lower 
branches, respectively. 

Compared with NR95, the most relevant study to ours, 
we have extended the expressions pertaining to slow 
waves in the slender slab limit (Eq.(16) in NR95) to ac- 
commodate negative values of phase speeds (see Eq. p^ V 
In addition, we have examined the case of wide slabs 
{kd ^ 1, see Eg. dTil) ) for slow waves. As for the fast 
waves, Eg. pT)) extends Eq.(13) in NR95 in the sense that 
only positive phase speeds are examined in NR95 but the 
phase speed in the present study is allowed to be nega- 
tive. Besides, we have presented the expressions of the 
dispersion behavior in the slender slab limit for the first 
fast kink branch (Eq.([T8])), and obtained new expressions 
for both kink and sausage waves at fed » 1 (Eg. lfTH)) ). 



3.2. Procedures for Computing Standing Modes 

To begin with, consider a pair of waves corresponding 
to {ijj,ki) and {uj,kr), where kr > ki in an algebraic sense. 
The combination of the two results in a Lagrangian dis- 
placement ^x{x^ z;t) in the form 

^x{x,z;t) = Re|f2;,i(a;)exp[-i(wt- fc(z)]| 

+ Re |^2:,r(a;) exp [—i {ujt — krz)]^ . (20) 

By "standing", we require that ^x{id,0;t) — and 
(,x{id, L;t) — 0, where z = and L denote the slab 
ends. Putting x — d and z = in Eg. ((20)) . one finds 

^x,i{d) — —^x,r{d), and hence it makes sense to choose 
£,x,i{d) — to be real. It then follows that 

^x{d,z;t) 

= [cos {ut — kiz) ~ cos (wi — krz)] 

. fkr-ki \ . ( ki + kr \ 
= —2A^ sm I z j sm lujt z j . (21) 

For ^xid, L] t) to be zero at arbitrary t, this requires 
27m 

fcr-fci = — ,n = l,2,--- (22) 

By convention, n = I corresponds to the fundamental 
mode, and n = 2 to its first overtone. 

Equation ()22|) suggests a graphical means to determine 
the period ratio P1/2P2 for a given aspect ratio d/L. 
This is illustrated by Figure [3] where we use the results 
computed for Mq = 0.8 as given in Figl^b. Here panels 
(a) and (b) are for kink and sausage modes, respectively. 
First, we convert the Vph — k diagram to an w — fc dia- 
gram. Only the first and second quadrants are relevant, 
with the former (latter) coming from the curves in the 
upper (lower) half-plane of Figl^b. Note that instead of 
combinations of positive k with negative Vph (and hence 
negative uj), we obtain combinations of negative k with 
positive Lu by reversing the sign of k in view of symmetry 
1. This, together with a negative Uph, yields a positive uj. 
We then compute P1/2P2 for a given d/L by measuring 
the width between the intersections of a horizontal line 
with the two curves corresponding to the waves that are 
to form standing modes. Take FiglSji for instance. Let A 
and B denote the intersections where a horizontal dashed 
line crosses bFKl and fFKl, respectively. The angular 
frequency wi for the fundamental mode is then found by 
requiring ABi = 2'Kd/L, while UJ2 for its first overtone is 
found by requiring that AB2 — A-Kd/ L. The period ratio 
is then P1/2P2 = uj2/2uji. In addition, from Fig|3)D one 
can readily see that the minimum allowed aspect ratio, 
namely {d/ L)cutof{, is not determined by the difference 
between the two cutoffs of the two uj — k curves involved 
divided by 27r, but larger than that. We will come back 
to this point in discussing FiglT] Furthermore, from sym- 
metry 2 one finds that if Mq is —0.8 instead, the uj — k 
diagram will be a refiection about the vertical axis of 
what we have for Afo = 0.8. It is then clear that the 
period ratios do not depend on the sign of Uo but only 
on its magnitude. 

At this point we note that it is impossible for a kink 
and a sausage wave to combine to form a standing mode. 
This is not obvious by examining Eg.p2p alone but has 
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to do with the symmetric properties of the eigenfunc- 
tions. Let kr and fc; represent, say, a sausage and kink 
wave, respectively. Now that the displacement for kink 
(sausage) modes is an odd (even) function of x, in light 
of Eq-dll]) one finds that 

^xi~d, z; t) — [cos {ujt — kiz) + cos {ujt — k^-z)] 

(kr-h \ ( ki+kr \ 

— zAf cos z cos ujt z , 

which does not have a permanent node at z = or L 
when ki and kr satisfy Eg. ip^ . hence violating the re- 
quirements for a standing mode. 

In computing the coronal standing modes, we consider 
only bFKl and fFKl for kink modes, and bFSl and fFSl 
for sausage modes. Higher order branches like bFK2 or 
bFS2 are not considered since they would form stand- 
ing modes only for relatively thick slabs. Take fFK2 and 
bFK2 for instance. For them to combine to form a stand- 
ing mode, the slab aspect ratio d/L has to exceed 0.85. 
In the case of fFSl and bFS2, d/L has to be larger than 
0.95. 

Can slow and fast kink (or sausage) waves form a 
standing mode? For sausage modes like fSS or bSS to 
combine with modes like fFSl or bFSl, the slab has to 
be very wide. In the case of fSS and bFSl, the aspect 
ratio d/L has to exceed 0.72. For kink waves, combi- 
nations like bFKl with fSK seem possible but turn out 
extremely unlikely since the slow mode tends to be dom- 
inated by the intensity oscillations rather than the dis- 
placement of the slab, whereas the opposite is true for 
the kink mode. For slender slabs kd — > 0, the backward 
fast kink wave corresponds to Vph ~ —vac, and nod — ^ 0, 

and hence X w , ^^^^\r^/i' — j(nQC?)^. If one neglects the 

flow and notes that c§ <C v\q, then one finds that this 
ratio is ~ {pq/ pe — l){kd)'^, meaning that fast kink waves 
are dominated by the displacement rather than the den- 
sity variation. However, for slow kink waves, by noting 
that for kd — > 0, (uph — Uq)"^ — > c^q, one finds that 
X « — (t;^Q/c§)(no(i) tan(nod). With nod approaching 
{j — l/2)7r, this will be large, meaning that slow kink 
waves are dominated by density variations. The end re- 
sult is that, if a fast kink wave does combine with a slow 
one to form a standing mode, a transverse displacement 
on the order of the slab width will cause a relative den- 
sity variation, and hence a relative intensity variation, 
that may exceed unity! This makes such combinations 
unrealistic. 

3.3. Period Ratios for Standing Kink Modes 

FigureSlpresents the period ratio P1/2P2 as a function 
of the slab aspect ratio d/L for standing fast kink modes. 
A series of values for the flow speed Uo is adopted, and 
Uq is measured in units of the internal Alfven speed 
Uq = Mavao- Common to all the curves is that they 
start with unity at small d/L, where wave dispersion is 
negligible, decrease with increasing d/L and attain a min- 
imum before rising again. With increasing Uq, the devia- 
tion of P1/2P2 from unity is enhanced substantially when 
compared with a static case Ma = 0. Take the minima 
for instance. While in the static case it reads 0.851, when 
Ma = 0.8 it is significantly reduced to 0.657, amounting 
to a relative difference of 22.8%. One may notice that 



the former is attained at d/L — 0.137, while the latter at 
d/L = 0.139. However, this is not to say that for smaller 
aspect ratios the effects of enhanced flow speed is negligi- 
ble. For instance, for d/L as small as 0.03, while P1/2P2 
reads 0.939 in the static case, it reads 0.85 when Ma in- 
crea ses to 0.8. From this we c onclu de that, in agreement 
with lMacnamara fc Roberts! ()2011[ ). transverse structur- 
ing may contribute to the deviation from unity of P1/2P2 
even for thin slabs. And we note that, in addition to the 
density contrast w hich is the sole contrib utor to wave 
dispersion in Macn amara fc Roberts! (j201lD , the velocity 
shear also plays a significant role. This has significant 
bearings on the applications of look-up diagrams express- 
ing the period ratio as a function of longit udinal density 
stratification like the one constructed by lAndries et al.l 
([2005) for coronal loops. While a slab diagram has yet to 
appear, it seems appropriate to say before using this di- 
agram to deduce the longitudinal stratification, one had 
better exclude the wave dispersion brought forth by the 
transverse structuring. This is especially necessary for 
strong transverse density contrasts, and/or in the pres- 
ence of substantial flow shears. 

The effect of flow is further shown in FiglSJ where (a) 
the minimal period ratio, (Pi/2P2)min, as well as (b) 
its location, (ci/L)inin, are shown as a function of the 
Alfvenic Mach number Ma- In addition to the nomi- 
nal value of VAe/vAo being 2, other ratios of 3 and 4 are 
also examined and shown in different colors. Consider 
Fig [5^ first. One can see that the effect of fiow on the 
period ratios is significant for all the considered VAe/vAo, 
or equivalently, the density contrast. Take vab/vao — 4 
for instance. While the effect of flow is slightly less pro- 
nounced than in the case where VAe/'^Ao = 2, introduc- 
ing a flow of Ma = 0.8 also rather significantly reduces 
the period ratio from the static value 0.778 to 0.633, or 
by 18.3%. Likewise, in the case with vab/vao = 3, a 
relative reduction of 20.2% is found with P1/2P2 being 
0.802 for Ma = but 0.64 for Ma = 0.8. Now consider 
Fig[5jD. One may notice that, with increasing density 
contrast, the aspect ratio at which the minimum period 
ratio is attained decreases. For a given density contrast, 
with increasing fiow this aspect ratio tends to increase, 
even though the effect is marginal for vac/vao = 2 and 
is more pronounced with higher vae/vao- For instance, 
for VAe/vAo — 2, a, value of 0.138 can be quoted for 
all the (d/L)inin obtained, whereas for vac/'^ao = 4, 
(<i/L)mi„ increases substantially from 0.057 in the static 
case to 0.106 when Ma ~ 0.8. Furthermore, the over- 
all tendency for both (Pi/2P2)min and (d/i)min to de- 
crease with i ncreasing density contrast ag rees closely 
with Fig. 10 in i Macnamara fc Robert j (|2011f ). However, 
while the minimum in the static case cannot go below 
\p2/2 — 0.707 as determined analytically for an Ep- 
st ein density profile and numeri cally for a step function 
in IMacnamara fc Robert j (|2011l ). it is not subject to this 
constraint in a flowing case. 

3.4. Period Ratios for Standing Sausage Modes 

Now consider standing sausage modes. Figure [6] 
presents the period ratio P1/2P2 as a function of aspect 
ratio d/L for both the cases where vac/vao = 2 (solid 
curves) and vac/vao = 4 (dashed curves). A series of 
Uq, once again expressed in terms of the Alfvenic Mach 
number Ma, is adopted and indicated by different colors. 
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For the case va^I'vaq = 2, one can see that except for the 
results corresponding to the first three values of Ma , the 
curves are absent since the range of d/L in which stand- 
ing sausage modes are allowed is outside the presented 
d/L range. For the three curves that are present, while 
one can hardly tell them apart when their overlapping 
parts are concerned, it is evident that they start with 
considerably different values of the allowed d/L. With 
Ma increasing from to 0.1 to 0.2, this cutoff aspect 
ratio {d/L)cutoS increases from 0.3 to 0.33 to 0.39. This 
means that at a given vac/vaoi relative to the static case, 
slabs with flow can support standing sausage modes only 
when they are sufficiently wide. The same is also true 
for other choices of vab/vao- when vac/vao = 4, while 
(d/L)cutoff reads 0.134 for the static case, it reads 0.157 
when Ma = 0.2 and 0.232 when Ma = 0.4. In this 
case, the effect of flow on the period ratio turns out to 
be more pronounced than in the case vac/vao = 2. At 
d/L = 0.35, the case with Ma = 0.5 yields a P1/2P2 of 

0. 654, which is 5.2% lower than the value 0.69 obtained 
in the static case. If comparing the solid and dashed 
curves, one can see that the period ratio as well as the 
cutoff aspect ra tio depend rather sensit i vely o n vac/vao, 
consistent with iMacnamara fc Roberts! (|2011l ). Here we 
have shown the new result that, for standing sausage 
modes, the flows are important in determining the ranges 
of aspect ratios where standing modes are allowed. For 
the considered density contrasts the relative changes in 
the period ratio due to flows are not as significant as in 
the standing kink cases, though. 

Before proceeding, we note in passing that due to the 
strong dispersion accompanying sausage waves, it is nat- 
ural to expect the perio d ratio Pi /2P2 to devia t e from 

1 . As demonstrated by IMacnamara fc: Roberts! (|2011l ) , 
for static coronal slabs, P1/2P2 may reach values as low 
as 1 /2 when the density contrast increases towards infin- 
ity. How low P1/2P2 is allowed to go for coronal loops 
with finite aspect ratios will be explored in a future study 
th at extends the sl ab computations presented here and 
in IMacnamara fc"^ obcrts (2011). Regardless, it is safe 
to say that P1/2P2 will deviate in general from 1 even in 
the absence of any longitudinal structuring since sausage 
waves are also strongly dispersive in a cylindrical ge- 
ometry. This implies that practices employing standing 
sausage modes to inf er longitudinal density scale heights 
may be problematic (Sriva stava et al.ll2008l ). 

Figure [7] further examines the effects of flow speed on 
the cutoff aspect ratio {d/L)cutoB pertinent to the stand- 
ing sausage modes. Here three cases with different values 
of VAe/vAo are examined, and presented with different 
colors. Besides, the dashed curves represent the distri- 
bution with Ma of the lower limit of (d/i)cutoff given 
by (A+ -t- A-)5i/27r = (A+ + A-)/4 (see Eq.l^). It 
is evident that with increasing Ma , the deviation of the 
computed values from this limit increases substantially. 
Note that because Cq « c^q ^ v\q, this limit can be 
satisfactorily approximated by 

which has only a modest Ma dependence in the parame- 
ter range considered. The fact that the computed cutoff 



aspect ratios have a much stronger Ma dependence re- 
sults from an increasingly strong asymmetry in the lo — k 
diagram, as discussed regarding FiglHJj. From an obser- 
vational viewpoint, the curves presented in Figl7]can be 
tested in a statistical sense if there is a large set of slab- 
supported standing sausage oscillations. Sorting the ob- 
servations into groups with similar density contrasts, and 
in any group constructing a parameter space comprising 
the aspect ratio as well as the Alfvenic Mach number, 
then one expects to see a ridge across which detections 
of such standing modes become drastically less frequent. 

A few words are in order to justify our examination 
of standing sausage modes for slabs with a finite width. 
Indeed, loops in the form of bright structures appear to 
be thin in most UltraViolet images of the solar corona, 
meaning that thin-tube limits are often valid in many 
seismological applications. However, coronal loops with 
finite aspect rati os may exist. For insta nce, the flaring 
loop examined in lNakariakov et aLl ()2003f) is 25 Mm long, 
and 6 Mm in diameter, resulting in an aspect ratio of 0.12 
using our definition. That loop was shown to be wide 
enough to support a global sausage mode. However, one 
may ask how typical these "fat" structures are, given that 
their observations seem to be rather sparse. With this 
caveat in mind, the results found here can nevertheless 
serve as a reference for our photospheric computations 
to compare with. 

4. PERIOD RATIOS FOR STANDING MODES SUPPORTED 
BY PHOTOSPHERIC SLABS 

4.1. Overview of Photospheric Slab Dispersion 
Diagrams 

Now move on to the photospheric case, by which we 
mean the ordering va o > Ce > cq > VAe h olds. Fol- 
lowing NR95 and Te rra-Homem et"al1 (|2003[ ). we deal 
only with an isolated slab immersed in an unmagne- 
tized atmosphere: cq = 1,wao = 1.5, Cg — 1.2, vac — 
(and hence cto = 0.83, ctc = 0,/9e/po — 2). Putting 
some absolute value for cq as 8 km s"^, one finds 
VAQ = 12, Ce = 9.6 km s^^, which agree with the photo- 
spheric values adopted by NR95, and fall in category (ii) 
in Evans & R oberts (1990) . 

Figure [5] presents, in a format similar to Fig HI the 
phase speeds Vp\i as a function of longitudinal wavenum- 
ber k for a series of Uq. Here again Uq is measured 
in units of the internal sound speed Mq = Uq/cq, and 
kink (sausage) modes are given by the dashed (solid) 
curves. Labeling different waves is slightly more com- 
plicated than in the coronal case, since now in addi- 
tion to the body waves, surface waves are also per- 
mitted and can be further grouped into fast and slow 
ones (Evans & Roberts 1990). However, given that body 
waves cannot combine with surface ones to form stand- 
ing modes (see section H?^ . and propagating body waves 
show only a modest dispersion, the flow effect turns out 
to be marginal for pure standing body modes. Therefore 
we leave them out, and label surface waves only, as indi- 
cated by FigOi. With the convention "b/f-hF/S-|-K/S" 
we have, for instance, fFS for forward fast sausage waves, 
and bSK for backward slow kink waves. 

Figure[5]indicates that the dispersion diagrams possess 
a stronger dependence on the internal flow than in the 
coronal case, evidenced by the fact that some propagat- 
ing windows may disappear. Take the body waves for 
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example, for which two bands exist and are shifted up- 
wards with increasing Uq ■ For forward (backward) ones, 
the phase speeds ?;ph start with cto + Uq {—cto + Uq) at 
fed ^ 1 in the same manner as described by Eq. ([T2l) , for 
both kink and sausage waves ahke. With increasing kd, 
Uph approaches ±co + Uo in the way given by Eg. (fT4l) . the 
only exception is the first branch of kink waves: Wph may 
exceed co + Uq (the first dashed curve from top in Fig|8^) 
or fall below — co + C/q (the first dashed curves from bot- 
tom in all panels), thereby rendering these body waves 
a surface one asymptotically. This phenomenon, termed 
mode crossing, was also shown in iTerra-Homem et alJ 
(|2003D . When - cro > Uq > Ce - cq (FiglHt), the 
upper window is bounded from above by Ce instead of 
Co -l- C/q, therefore providing a maximum allowed longi- 
tudinal wavenumber beyond which forward slow body 
kink waves are no longer trapped. When Uq > Ce — cto 
(FiglHJi), this upper window disappears altogether. 

Surface waves are examined hereafter. Consider first 
the fast ones. One finds that for kd <^ 1, only sausage 
waves are allowed, kink waves appear only at relatively 
large kd as a continuation of the body kink waves, i.e., 
via mode crossing. This is readily understandable given 
that with kd ^ and v^^^ — >• c^, if mod then the DR 
for kink waves (Eq.©) does not allow any solution since 
its LHS approaches infinity. For the fast sausage wave, 
one has 
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where Ce.± = ±Ce — Uq. Now consider the slow surface 
waves. One can readily see that for sausage waves Uph 
starts with ±cto + Uq, and for small fc. 



Vph ~ ±cto 
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where Ctq.± = ±cto + C^- At sufficiently strong internal 
flow with Uo > Ce — CTo, the upper branch of slow sausage 
waves starts with Ce instead, and the behavior of Wph at 
small k is identical to Eq. (P5)) with the plus sign. On 
the other hand, for kink waves, Wph at small k can be 
approximated by 



Vph 



± 



PO , 2 



AO 



U'^)kd. 



(25) 



The study so far conducted in this section differs sub- 
stantially from NR95 (their section 3.3), even though 
both studies are for isolated photospheric slabs. NR95 
examined the consequences of a variable Ue while keeping 
= 0, as opposed to our study where t/g = but Uq is 
allowed to vary. A situation similar t o our choice was ex- 
amined in lTerra-Homem et al.l (|2003f ) albeit a cylindrical 
geometry is considered there. 

4.2. Computing Standing Modes 

As have been discussed, when constructing standing 
modes, we do not employ pure body waves in view of 
their modest dispersion. Once again, by saying it is pos- 
sible or impossible for a pair of propagating waves to form 



a standing mode, we mean the standing mode hence con- 
structed corresponds to a realistic density fluctuation X 
in the slender slab limit kd ^ provided that the waves 
have no cutoff longitudinal wavenumber. 

Consider kink waves for kd ^ first. To this end, 
let us recall that fast kink surface waves do not ex- 
ist, whereas for slow kink surface waves one finds that 
X « [/q/(c§ — UQ){mod)'^ by noting that Wph — >■ and 
mod — > 0. On the other hand, the kink body waves in 
the limit fed —> is proportional to (nod) tan(no(i) which 
approaches infinity given that ngd — >■ (j — l/2)7r. This 
means that standing kink modes can result from a com- 
bination of either pure body or pure surface waves, but 
not a mixture of both. Therefore we are left with only 
one task of examining pure kink surface modes. 

Standing sausage modes can be considered in a similar 
manner. Once again, they cannot be formed by combin- 
ing a pair of body and surface waves. This is because 
for body waves, X is proportional to {ngd) cot{nod) and 
approaches infinity when kd approaches zero given that 
nod — jir. For surface waves, however, one finds that 
at kd ^ 1, for fast waves X ^ cI/{cq — c1), whereas for 
slow waves 



X 



if Uo < Ce - CTO 



(c^-Uof 



. if t/o > Ce 



(26) 



Cto 



Note that the latter holds since Wph — >■ Cg when Uq > Ce — 
Cto- Nonetheless, different from the coronal case, now in 
the surface category fast and slow waves are allowed to 
form standing sausage modes. 

4.3. Period Ratios for Standing Kink Modes 

Figure [9] presents the period ratio Pi /2P2 as a func- 
tion of the slab aspect ratio d/L for standing kink sur- 
face modes. A series of values for the flow speed Uo is 
adopted and given by different colors. Common to all 
curves is that they start with ~ 1.4 at small d/L and 
decrease monotonically with d/L towards unity. The 
behavior for slender slabs is in stark contrast to the 
coronal case where for extremely thin slabs d/L ^ 
the period ratio approaches unity no matter how strong 
Uq is. What happens in the photospheric case is, irre- 
spective of Uo, P1/2P2 differs substantially from unity 
even for zero- width slabs! This is understandable given 
that for small kd the surface kink waves have phase 
speeds Vph oc Vkd, leading to that the angular frequency 
uj — Vphk is proportional to k^^"^. Note that, as il- 
lustrated by FiglH the period ratios are computed via 
P1/2P2 — UJ2/2LO1, where uji and uj2 correspond to the 
angular frequencies found at the separation AB being 
2TTd/ L and twice that, respectively. This means, at small 
d/L, P1/2P2 « 2^/2/2 = V2 = 1.414. 

One can see that the effect of flow on P1/2P2 is sub- 
stantial for d/L only in the range between 0.03 and 0.27. 
Take d/L = 0.1 for example. While P1/2P2 = 1.24 in 
the static case, it reads 1.1 when Mq = 0.4, amounting 
to a reduction of 11.6%. When d/L < 0.03, the /c^/^ de- 
pendence of phase speed dominates, leading to little dif- 
ference between different curves. To the contrary, when 
d/L > 0.27, all curves show little deviation from unity 
since the wave dispersion is only modest (see Fig|8]). 



Period ratio of magnetic slabs 



9 



4.4. Period Ratios for Standing Sausage Modes 

Consider now standing sausage surface modes. Fig- 
ure [10] presents the period ratio P1/2P2 as a function of 
aspect ratio d/L for a series of Uq. Distinct from the 
coronal case, now there is no cutoff in d/L below which 
no standing modes can be constructed. In addition, the 
combinations of fast and slow waves are now allowed, re- 
sulting in four pairs of possibilities. Note that instead 
of four, there are only two curves in FigslTUb and [TUti. 
since forward fast waves are absent in the rest of the 
cases considered. In the parameter range we explored, 
the strongest Uq dependence of P1/2P2 happens in the 
case where backward fast waves combine with forward 
slow waves, presented in FigfTUb. The trend of P1/2P2 
with changing d/L is the same for all curves in this case, 
decreasing from unity at small d/L to some minimum 
and then increasing toward unity again. Even in this 
case, the flow effect is only modest, as evidenced by the 
fact that for Mq = 0.4, the minimum P1/2P2 attains is 
0.909, only slightly smaller than 0.965 attained in the 
static case. 

Directly measuring sausage oscillations with imag- 
ing in struments was no t possible until very re- 
centlv (iMorton et all 1 20111) (hereafter M EJM. also see 
iDorotovic et al.ll2008l iMorton et al.ll2012[) . By using the 



Rapid Oscillations in the Solar Atmosphere (ROSA) in- 
strument situated at the Dunn Solar Telescope, MEJM 
identified sausage modes in magnetic pores by examin- 
ing the phase relations between the pore size and its 
intensity. A number of periods were present, and the 
authors suggested that some of the values may be at- 
tributed to the overtones of the fundamental mode, with 
the standing mode set up by reflections between the pho- 
tosphere and transition region. Let Pn iji = 2,3, ..) de- 
note the period of the overtones with Pi standing for the 
fundamental. For now, let us consider only the period 
281 ± 18 s found in the intrinsic mode function (IMF) 
c4 therein. Choosing Pi to be 550 s pertinent to fast 
modes, this value would correspond to a Pi/nP„ in the 
range [0.92, 1.05] ([0.61, 0.70]) provided it corresponds to 
n = 2 (n = 3). Note that the aspect ratio of the struc- 
ture in question is d/L ss 0.3, since the average pore 
size is 2d « 1.16 x 10^ km (corresponding to an area 
of 1.36 X 10^ km2) and L = 2 x 10^ km. At this as- 
pect ratio, the computed P1/2P2 as given in our FigfTOl 
indeed lies in the range that corresponds to n = 2, re- 
inforcing the suggestion by MEJM that it is the first 
overtone. Nevertheless, choosing Pi to be 660 s perti- 
nent to slow modes, this period would correspond to a 
Pi/nPn in the range [1.1,1.26] ([0.74,0.84]) provided it 
corresponds to n = 2 (n = 3). While the former seems 
inappropriate, the latter is plausible given that the dis- 
persion introduced by transverse structuring in density 
and magnetic field strength, along with the possible role 
of velocity shear, reduces Pi/nPn from unity for pho- 
tospheric sausage modes. While this comparison is ad- 
mittedly inconclusive, the point we want to make is that 
the transverse structuring does broaden the range of the 
period ratios, thereby offering more possibilities in inter- 
preting the measured oscillation periods. Instead of pro- 
ceeding further with this discussion, let us note that the 
slab model may be seen as only a first approximation of 
solar atmospheric structures, and the oscillations mea- 



sured by ROSA may correspond to propagating rather 
than standing waves, as was pointed out by MEJM. 

5. SUMMARY AND CONCLUDING REMARKS 

The present study has been motivated by two 
strands of research in the field of solar magneto- 
seismology. On the one hand, exploiting the mea- 
sured multiple periods is playing an increasingly im- 
portant role in deducing the lon gitudinal structuring 
of solar atmospheric stru ctures (|Andries et al.l 120091 : 
iRuderman fc Erdelvil 120091 ) On the other hand, it was 
recently realized that significant flows in coronal struc- 
tures reaching the Alfvenic regime can bring signif- 
icant revisions to the seismologically deduced physi- 
cal parar neters such as the co ronal loop magnetic field 
strength (jTerradas et al.|[2011f) . To complement and ex- 
pand the only theoretical study on multiple periodicities 
in coronal slabs (Macnaniara & Roberts 2011), we ex- 
amine the consequences the flows may have. To this 
end, we numerically solve the dispersion relations for 
waves supported by slabs incorporating flows, and de- 
vise a graphical method to construct standing kink and 
sausage modes. 

For coronal slabs, we find that the internal flow has 
significant effects on the standing modes, for the fast 
kink and sausage ones alike. For the kink ones, they may 
reduce the period ratio by up to 23% compared with the 
static case, and the minimum allowed period ratio may 
fall below the analytically derived lower limit in the static 
case. In particular, the reduction due to a finite fiow may 
be significant even for thin slabs. For the sausage modes, 
while introducing the flow leads to a reduction in the 
period ratio that is typically < 5% relative to the static 
case, it has significant effects on the threshold aspect 
ratio only above which standing sausage modes can be 
supported. At a given density contrast, this threshold 
may exceed its static counterpart by several times for the 
parameters considered, thereby limiting the detectability 
of standing sausage modes to even wider slabs. 

For the isolated photospheric slabs, we find that the 
internal flow has only a marginal effect on the period 
ratios P1/2P2 for the surface modes, and even less for 
body modes. That said, standing modes supported in 
this case are distinct from the coronal case in a number 
of aspects. First, standing sausage modes may be sup- 
ported by slabs with arbitrary aspect ratios and do not 
suffer from a cutoff d/L any longer. Second, for standing 
kink surface modes, P1/2P2 deviates from unity even for 
a zero-width slab, which originates from the fc^^^ depen- 
dence of phase speeds of slow kink surface waves in the 
slender slab limit. 

The present study can be extended in several ways. 
First, to be able to directly compare with coronal loop 
oscillations, so far the most observed and the best un- 
derstood, we need to turn to the cylindrical geometry 
and examine how introducing the flow affects the disper- 
sion behavior and consequently period ratios and cutoff 
aspect ratios. Second, allowing the tube parameters to 
be time-dependent has been found important as far as 
the seis mological applications of t he period ratio are con- 
cerned (|Morton fc Erdelvill2009af) . This is expected to be 
also true when the flow speeds are dynamic rather than 
time-stationary. Third, concerning the period ratios, one 
important contributor to its deviation from unity that 
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has received much atten tion is the longitudinal structur- 
ing (|Andries et aL|[2009l ). which has yet to be considered 
in examining the slab period ratios. Nonetheless, the 
main results presented here, which demonstrated the in- 
fluences of the introduced flow on the standing kink and 
sausage modes, point to the need for further investiga- 
tions into such influences. 

Before closing, let us recall that solar magneto- 
seismology (SMS) has extended well beyond the terri- 
tory originally taken up by "coronal seismology" . One 
thing to note in particular, as advocated by Morto n~et al.l 
(|201lD . is to seriously consider applying SMS to photo- 
spheric structures because there exist ample observations 
of oscillatory motions on the one hand, relevant theories 
are available on the other. Strictly speaking, our results 
on the sausage modes are more applicable to the lower 
parts of the solar atmosphere. In this sense, our com- 



putations strengthen the notion raised by iMorton et al.l 
(|201lD that photospheric studies in the context of SMS 
offer equally rich physics as their coronal counterparts. 
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per substantially. This research is supported by the Na- 
tional Natural Science Foundation of China (40904047, 
41174154, and 41274176), the Ministry of Education of 
China (20110131110058 and NCET-11-0305), and by the 
Provincial Natural Science Foundation of Shandong via 
Grant JQ201212. 
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APPENDIX 

ISSUES RELATED TO SYMMETRY PROPERTIES OF THE DISPERSION DIAGRAMS 

In this appendix, we examine the two symmetry properties presented in section [5] in some detail. In particular, we 
show that for the purpose of examining the period ratios of standing modes supported by a flowing structure embedded 
in a static background, one needs to consider only positive longitudinal wavenumbers k and positive internal flow speed 

Figure [TT] illustrates what happens in a general case where neither Uq nor Ue is zero. It presents the wave phase 
speed Wph as a function of longitudinal wavenumber fc, in a format similar to Figsl2]and[8l Note that both positive and 
negative k-s are considered here. Kink and sausage waves are given by the dashed and solid curves. The horizontal 
and vertical dash-dotted lines in purple separate the graph into four quadrants. The parameters used to construct 
this graph are identical to those adopted in Figl^b, except that now the external medium is not static, Ue = cq to 
be specific. This rather arbitrary choice of Ue is to illustrate symmetry property 1, which in fact holds for arbitrary 
values of Ue and Uq. Graphically, this property translates into the symmetry between the left half-plane (the 2nd and 
3rd quadrants) and the right one (the 1st and 4th quadrants). In this sense one needs only to consider the 1st and 
4th quadrants as long as one is interested only in trapped waves. 

The second symmetry property is best illustrated if one compares Fig |12l with Fig ll31 where Vph as a function of k 
is shown for both a positive internal fiow (panel a) and a negative one (panel b) with |J7o| = O.Scq. The characteristic 
speeds are identical to those adopted when constructing Figl^l The difference between FigsH^] and [T3] is that in the 
former Ue is not zero (Ue chosen to be cq for illustrative purpose), whereas in the latter Ue = 0. To label the waves, 
the same convention as in Figsl2] and |8] is adopted. Note that in this coronal case no surface modes are introduced 
by the non-zero external flow, regardless of the sign of the internal fiow. The body waves can, as in the case where 
Ue = 0, be grouped into slow and fast ones. The phase speed for the former is bounded by ±cto + Uq and ±co -f Uq, 
and Vph for the latter is bounded by ±VAe + Ue and ±vao + Uq. The most important difference, as far as the symmetry 
properties are concerned, is that Fig |12l lacks a symmetry between the two panels, whereas in Fig |13l the curves in the 
1st (4th) quadrant of panel (a) map to those in the 4th (1st) quadrant of panel (b). For instance, the curve labeled 
bFKl in FigHSli corresponds to fFKl in FigHSb, and likewise fFK2 i n FigjlSh corresponds to b FK2 in FigUSt). As 
a matter of fact, this symmetry was present already in Figs. 4 and 5 in iTerra-Homem et al.l (|2003[ ) where wave modes 
supported by a fiowing cylinder is examined. This is not surprising given that the authors also examined the case 
where the external background is static, and that the Uq they adopted in constructing Fig. 4 {Uq = —0.25vao) is close 
to the one used for constructing Fig. 5 (Uq = 0.35?;^o)- A perfect symmetry would result if the authors chose a pair of 
J7o of the same magnitude but with opposite signs. 

The second symmetry, which exists only in the absence of external flows, docs not invalidate the concept of reso- 
nant fiow instalDility (RFI) (Rvutova 1988; HoUwcg ct al. 1990; Erdelyi & Goosscns 1996) (see Goosscns ct al. 20111 
and iTarovan fc RudermanI I2011i for recent reviews). This point can be illustrated if one examines in detail, say, 
lAndries et al.l (|2000f ) where the authors address the possible instability a velocity shear between a coronal plume and 
its environment may introduce. The authors choose to work in the zero-beta limit, and in a frame in which Ue = and 
Uq is negative. When a continuous Alfven speed profile connects the internal and external values, RFI sets in when 
the "originally forward-propagating" (in the sense that Wph is positive when Uq = 0) kink waves reach the continuum 
of the originally backward-propagating Alfven waves. This happens at a Uq smaller in magnitude than required by 
the non-resonant Kelvin-Helmholtz instability, the difference between the two thresholds being quite significant for 
high density contrasts between th e internal and external media. In light of the symmetry property we discussed, 
what Fig. 3 in lAndries et al.l (|2000[ ) shows for a negative Uq can help one to foresee what happens when Uq is positive. 
In that case, RFI will also happen when Uq reaches exactly the same magnitude, for at this point the originally 
backward-propagating kink waves reach the continuum of the originally forward-propagating Alfven waves. Extending 
the present study to address RFI is physically very appealing, but may require quite some effort given the complexity 
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related to solving eigenvalue problems that involve dealing with the dissipative layers at the Alfven resonances and 
the possible cusp resonances. 

Figure [14] shows the wave angular frequency w as a function of fc, constructed from the Vph — k diagram presented in 
Fig[T3l Here FigslTik and[T4b correspond to Uq = O.Scq and — O.Scq, respectively. Only kink waves are presented, for 
otherwise the graph will be too crowded. The first thing to note is that if one insists on cj being positive, one would 
need the curves in the 3rd quadrant of a Vph — k diagram where both k and Uph are negative. However, from symmetry 
property 1 illustrated in FigHTl one can construct the 3rd quadrant curves using the curves in the 4th quadrant by 
keeping Vph unchanged while reversing the sign of k. Second, Fig ll4l inherits the same symmetry from Fig ll3l For 
instance, the branch labeled bFKl in FigfMk maps to the one labeled fFKl in FigfT4b. Likewise, fFKl in Fig|T4k 
corresponds to bFKl in FigfMb. Recall that to compute the period ratios for a given aspect ratio d/L, we simply 
choose the two curves corresponding to the propagating waves that are to combine to form a standing mode, and 
then measure wi and uj2 that correspond to the two horizontal cuts resulting in a separation being 2TTd/L and twice 
that, respectively. It follows from this procedure that the period ratio thus calculated would be the same in FigsfHk 
and 114b . Evidently, if we choose to work in a frame where the external medium is not static, then the symmetry 
between FigsfT4k and [T4b will be broken, meaning that the period ratio for the same d/L derived for a positive Uq 
will differ from that derived for a negative one, even if the magnitudes of Uq are equal. This means, at any given 
non-zero Ue, for one to examine the effects of Uq on the period ratios one would have to first calculate a full range 
of different Uq^ both positive and negative, then construct the uj — k diagram, and finally measure wi and UJ2- This 
certainly merits a detailed discussion, but is beyond the scope of the present study. 
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Fig. 1. — Schematic diagram illustrating the magnetic slab (denoted by subscript 0) and its environment (subscript e). The variables 
Pi,Ci,VAi and Ui {i = 0, e) represent the mass density, adiabatic sound speed, Alfven speed, and the flow speed, respectively. 
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Fig. 2. — Phase speeds Dpjj as a function of longitudinal wavenumber k for a series of internal flow Uq. Expressing Uo in units of the 
internal sound speed Uq = Mqcq, panels (a) to (d) correspond to an Mq of 0,0.8,1.2 and 3.2, respectively. On the left (right) of each 
panel, the characteristic speeds external (interior) to the slab are given by the horizontal bars. Kink and sausage modes are presented by 
the dashed and solid curves, respectively. They are further labeled, as showrn in panel (b), using combinations of letters b/f+F/S+K/S, 
representing Backward or Forward, Fast or Slow, Kink or Sausage. The number appended to the letters denote the order of occurrence. 
Hence, bFKl represents the first branch of backward Fast Kink mode. Moreover, here = 8, = 0.72, cj-^ = 0.719, while Vj^q = 4, 
Co = 1, and cxo = 0.97. 
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Fig. 3. — Determination of the period ratio with ui — k dispersion diagrams for (a) standing kink and (b) standing sausage modes. Different 
branches labeled with combinations of letters and numbers follow the same convention as in Fig[2] 




Fig. 4. — Period ratio P1/2P2 as a function of the slab aspect ratio djL for standing fast kink modes. Curves with different colors 
represent results computed for different values of the flow speed C/q measured in units of the internal Alfven speed t/o = M^v^a. 
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Fig. 5. — Effects of flow speed on (a) the minimal period ratio, (Pi/2P2)mini and (b) the aspect ratio at which the minimum is attained, 
(d/L)min. Here both (Pi/2P2)min and (d/-L)min are displayed as a function of the Alfvenic Mach number Ma- 



Period ratio of magnetic slabs 



17 





n 


o u 







75 







70 














\ 





65 










U 


a Pi 
oU 







55 







50 



0.1 





0.1 

0.2 
0.3 
0.4 
0.5 



0.2 



0.3 

d/L 



0.4 



0.5 



Fig. 6. — Period ratio P1/2P2 as a function of the slab aspect ratio d/L for standing sausage modes. The solid and dashed curves are for 
the cases where VAe/^AO = 2 and 4, respectively. Curves with different colors represent results computed for different values of the flow 
speed Uo measured in units of the internal Alfven speed Uo = Mj^Vj^q. 




Fig. 7. — Effects of flow speed on the cutoff aspect ratio, (a!/ij)cutoff > the lowest aUowed for standing sausage modes to occur. Here Mji 
measures tlie internal flow speed Uo in units of the internal Alfven speed, and different colors are for different ratios of VAe/^AOy ranging 
from 2 to 4. The dashed curves are the lower limit for (rf/ij)cutoffi accurate for the static case (please see text for details). 
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Fig. 8. — Similar to Fig|2]but for an isolated photospheric slab (vAe = 0). The characteristic speeds are Ce = 1.2cq,vao = I-Scq and 
cto = 0.83co. Panels (a) to (d) correspond to an Mq of 0,0.1,0.25 and 0.4, respectively, where Mq = Uq/cq measures the internal flow in 
units of the internal sound speed. Body waves, namely those with phase speeds between ztc^o + Uq and ±co + C/o, are not labeled. Surface 
waves are labeled with combinations of letters b/f+F/S+K/S, representing backward or forward. Fast or Slow, Kink or Sausage. 




Fig. 9. — Period ratio P\l2Pi as a function of the slab aspect ratio djL for standing slow surface kink modes. Curves with different 
colors represent results computed for different values of the flow speed C/q measured in units of the internal sound speed (7o = Afoco- 
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Fig. 10. — Period ratio P1/2P2 as a function of the slab aspect ratio d/L for standing surface sausage modes. Presented are combinations 
of (a) backward slow + forward slow, (b) backward slow + forward fast, (c) backward fast + forward slow, and (d) backward fast + 
forward fast waves. Curves with different colors represent results computed for different values of the flow speed Uo measured in units of 
the internal sound speed Uq = Mqcq. 
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Fig. 11. — Phase speeds tiph as a function of longitudinal wavenumber fc in a general case where neither Uo nor Ue is zero. The 
characteristic speeds are the same as adopted in FigjS] namely, = 8, Ce = 0.72, Vj^q = 4, and cq = 1. For Ue and Uo, we choose them 
to be Co and O.Scq, respectively. On the left (right) of the figure, the Doppler-shifted characteristic speeds external (interior) to the slab 
are given by the horizontal bars. Kink and sausage modes arc presented by the dashed and solid curves, respectively. The horizontal and 
vertical dash-dotted lines in purple split the graph into four quadrants. Note the symmetry between the left and right halves of the graph 
with respect to the vertical purple line, i.e., the fph— axis. Actually this symmetry is valid for arbitrary Uo and Ue. 
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Fig. 12. — Phase speeds tiph as a function of longitudinal wavenumber k where Ue is non-zero. The characteristic speeds are the same as 
adopted in Fig[2l namely, v^^ = 8, Ce = 0.72, v^q = 4, and cq = 1. Panels (a) and (b) correspond to a Uo of O.Scq and — O.Scq, respectively. 
In both panels, Ue = cq. On the left (right) of the figure, the Doppler-shifted characteristic speeds external (interior) to the slab are given 
by the horizontal bars. Kink and sausage modes are presented by the dashed and solid curves, respectively. The naming convention for 
the wave modes is the same as in Fig[2] 
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Fig. 13.— Similar to Fig[l2]but for the case where Ue = 0. Note the symmetry that maps a curve in the 1st quadrant in (a), the one 
labeled fFKl for instance, to a curve in the 4th quadrant in (b), the one labeled bFKl in this example. Likewise, this symmetry maps the 
curves in the 4th quadrant in (a) to the ones in the 1st quadrant in (b). This symmetry takes place only when Ue is zero. 




Fig. 14. — Angular frequency as a function of k for propagating kink waves supported by a slab embedded in a static medium {Ue = 0). 
Pa nels (a) and (b) correspond to a Uq of O.Scq and — O.Scq, respectively. Note the symmetry between (a) and (b), which is inherited from 
Fig |13l Consequently, this symmetry happens only when Ue = 0. 



